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MnSi is a cubic compound with small magnetic anisotropy, which stabilizes a helimagnetic spin 
spiral that reduces to a ferromagnetic and antiferromagnetic state in the long- and short-wavelength 
limit, respectively. We report a comprehensive inelastic neutron scattering study of the collective 
magnetic excitations in the helimagnetic state of MnSi. In our study we observe a rich variety 
of seemingly anomalous excitation spectra, as measured in well over twenty different locations in 
reciprocal space. Using a model based on only three parameters, namely the measured pitch of the 
helix, the measured ferromagnetic spin wave stiffness and the amplitude of the signal, as the only 
free variable, we can simultaneously account for all of the measured spectra in excellent quantitative 
agreement with experiment. Our study identifies the formation of intense, strongly coupled bands 
of helimagnons as a universal characteristic of systems with weak chiral interactions. 



PACS numbers: 

I. INTRODUCTION 



The spontaneous breaking of a continuous symmetry 
in a magnetically ordered state implies the existence of 
Goldstone modes. Thus, the low-energy spin wave dis- 
persions of ferromagnets and commensurate antiferro- 
magnets are universal - they can be deduced from simple 
symmetry arguments and the description does not de- 
pend on microscopic details 1 -"^. In ferromagnets, where 
the order parameter is a conserved quantity, the resulting 
spin waves have a quadratic dispersion while in antifer- 
romagnets the superposition of the normal modes of the 
sublattices leads to the well-known linear dispersion at 
low energies^'. 

In recent years complex forms of magnetic order re- 
lated to weak chiral interactions of non-centrosymmetric 
systems, also referred to as Dzyaloshinsky-Moriya 
(DM) interaction attract great interest. Systems in 
which DM interactions play an important role in- 
clude multiferroics^, the parent compounds of the high- 
temperature superconductors™, thin magne tic fil ms (Mn 
on a W substrate, or Fe on a Ir substrat e! 10 * 11 !), heavy 
fermion superconductor^^ and even itinerant-electron 
magnets, like MnSi, which displays a skyrmion lattice 
at ambient pressure^ and an extended non- Fermi liquid 



resistivity at high pressure d 14 * 15 !. As DM interactions are 
of growing importance in non-centrosymmetric systems 
this raises the question for any universal properties of the 
spin excitations in systems with weak chiral interactions. 
Helimagnetic order is very well established in systems 
with magnetic frustration such as rare earth elements like 
Tb or Ho. An excellent experimental and theo retica l ac- 
count has been given by Jensen and Mackintos h^ 6 ! 17 !, but 
excitation spectra perpendicular to the magnetic propa- 
gation vector have not been studied before. In addition it 
has long been appreciated that helimagnetic order due to 
frustration differs fundamentally from helimagnetic or- 
der due to DM interactions.^ The work on frustrated 
magnets has been contrasted by studies of the tetragonal 
system Ba2CuGe207E2l in which the influence of DM in- 
teractions on the excitation spectra has been addressed. 
However, Ba2CuGe20y orders antiferromagnetically on 
local length scales. Due to the much stiffer spectrum of 
an antiferromagnet compared to a ferromagnet, the ef- 
fects of DM interactions in a helical antiferromagnet are 
therefore substantially weaker and qualitatively different 
as compared to a helical ferromagnet like MnSi. More- 
over, the helimagnetic order in Ba2CuGe2C>7 is distorted 
and shows higher order harmonic contributions of at least 
20%. 

We have therefore decided to study the perhaps simplest 
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FIG. 1: (a) Depiction of a low-energy spin wave excitation of 
a ferromagnet with wavelength A cx . The spins precess around 
the uniform magnetization, (b+c) Sketch of a similar excita- 
tion in the case of a helimagnet with a momentum, q, parallel 
to the propagation vector of the helix, kh. The wavelength 
A cx is smaller (b) or much larger (c) than the pitch of the 
helix Ah = Note that only for this parallel configuration 
the excitations of the ferromagnet and the helimagnet are 
very similar for short A cx . (d) Reciprocal space map around 
the nuclear (110) Bragg peak. The locations of the magnetic 
Bragg satellites ki to k4 associated with the four domain 
configurations of the helimagnetic order are shown as color 
coded throughout the text. The locations at which the spin 
excitations were measured are shown by the open symbols. 
They may be grouped along three trajectories, where the res- 
olution ellipsoids are shown in gray shading. Note that all of 
the spectra recorded may be simultaneously accounted for by 
the model described in the text using just one parameter (the 
intensity) . 



example for the effects of DM interactions, notably cu- 
bic chiral helimagnets, where the uniform magnetization 
rotates slowly around an axis with a characteristic wave 
vector k/j. As illustrated in Fig.[TJ a helical state is es- 
sentially a ferromagnet on short length scales and an an- 
tiferromagnet for long distances. In fact, from a more 
general viewpoint, all forms of complex order can be in- 
terpreted as a superposition of helimagnetic order^Sl. In 
cases where the pitch of the helix is much larger than the 



lattice spacing, A/j ^> a, a helix can become remarkably 
stable against crystalline anisotropies. The Goldstone 
theorem ensures that the excitations of an incommen- 
surate helimagnet are gapless (in case the heli magn et is 
commensurate the gap is exponentially small) ' '. In 
this paper we report, that the spin excitations for a he- 
lix with a long pitch display an universal property: the 
formation of strongly coupled bands of helimagnons. 
The B20 compound MnSi is ideally suited to study the 
collective spin excitations of helimagnets due to DM in- 
teractions experimentally, because the magnetic prop- 
erties result from a clear separation of energy scales 
in a metallic host. MnSi crystallizes in the non- 
centrosymmetric cubic space group P2i3 (a = 4.558 A). 
Below T c = 29.5 K and in zero magnetic field a long- 
wavelength spin spiral with the spins perpendicular to 
the propagation direction stabilizes. The periodicity 
A/j ~ 180 A of the helix results from the competition 
of ferromagnetic exchange interactions, as the strongest 
scale, and Dzyaloshinskii-Moriya (DM) interactions as 
a manifestation of weak spin- orbit coupling in crystal 
stru ctures without inversion center, on an intermediate 
scal c! 23 * 2 ^ . The propagation direction of the spin spi- 
ral is finally locked to the cubic space diagonal through 
very small crystal field interactions, providing the weak- 
est scale. In comparison with the helical modulation, the 

Fermi wave- vector is large kp ~ 0.7 A 
MnSi has recently attracted great interest as the per- 
haps best candidate displaying a genuine non-Fermi liq- 
uid met allic sta te in a three-dimensional metal at high 
pressure^ 14 * 15 * 25 !. Moreover, based on neutron scattering 
and /.i-SR at high pressure s ' 26 1 27 1 , it has been suggested 
that the NFL behavior may be related to spin textures 
with non-trivial topologjEHEil j n fact, a skyrmion lat- 
tice, was recently identified unambiguously at ambient 
pressure in a small phase pocket just below T c , which is 
believed to be stabilized by thermal fluctuations^. To 
resolve the origin of these exciting properties it has be- 
come of great importance in its own right to establish at 
first a full account of the normal helimagnetic state and 
its spectrum of excitations. This has also inspired the 
experimental study presented here. 

The three different energy scales governing the physics 
of the helimagnetic state of MnSi result in four different 
regimes for spin excitations. As a function of increas- 
ing momentum of the excitations these regimes may be 
summarized as follows. For the smallest momenta one 
explores the antifcrromagnetic limit, q <C k p , where a 
linear spectrum of spin waves is expected. Here q is 
the wavevector measured from the ordering vector kh of 
the helix and k p measures the weak pinning of the ori- 
entation of the helix to the crystalline lattice, which is 
fourth order in spin orbit coupling (k p is expected to be 
of the order of k\jkp and therefore extremely tiny). The 
antiferromagnetic limit is followed by the helimagnetic 
regime for k p <C q <C fch, where the pinning of the he- 
lix can be neglected. Here the universal properties are 
predicted to consist of a dispersion similar to an antifer- 
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romagnet, u> oc |q|" with n = 1 for propagation parallel 
to kh, whereas in the transverse direction the spectr um is 
quadratic with n = 2 characteristic of a ferromagnet P 1 * 3 2 1 
Further increasing the momenta leads to the cross-over 
between the helimagnetic regime and ferromagnetic limit, 
^ Q <C k c with k c ~ kp deep in the magnetic phase. 
The magnetic structure is now probed on a length scale 
short compared to its pitch and a spectrum reminiscent 
of a ferromagnet may be expected, because the system 
locally appears to be a ferromagnetically aligned. How- 
ever, as shown below this expectation differs dramatically 
from the observed behavior. Instead of ferromagnetic 
magnons multiple hclimagnon bands are excited, which 
may be regarded as a new universal property of magnetic 
materials. Finally, for the largest momenta, g > fc c , the 
ferromagnetic spin waves cross over into the Stoner con- 
tinuum as established experimentally in a large number 
of systems. The properties depend now on details of the 
Fermi surface, i.e., they are no longer universal. 
Several pioneering inelastic neutron scattering studies 
have been carried out in MnSi. In the mid 1970ies 
Ishikawa and coworkers established for the first time in 
any magnetic metal the existence of paramagnon fluctu- 
ations over major portions of the Brillouin zone 3 - 3 -'. They 
further observed ferromagnetic spin waves when sup- 
pressing the helimagnetic state in an applied magnetic 
fielcpS. More recently polarized inelastic neutron scat- 
tering established that the spin fluctuations in the para- 
magnetic state are chiraP^ and t he exi stence of spin-flip 
excitations in a small applied fiel d 35 ' 36 '. 
The magnetic phase diagram of MnSi is remarkably rich. 
For small magnetic fields between of 0.05 T and 0.15 T 
a reorientation transition takes place and the helimag- 
netic modulation aligns with the applied field to form the 
so-called conical phase. The helical modulation is sup- 
pressed for fields above ~ 0.6 T. Finally, a small phase 
pocket just below T c has been identified recently as the 
first example of a skyrmion lattice, thereby demonstrat- 
ing the inherent instability of the helical order to stabilize 
non-trivial spin textures in small magnetic fields. 
For a comprehensive understanding it is consequently im- 
portant to measure the spin excitations in the prestine 
(zero-field cooled) helical state. This has not been at- 
tempted before, since four equally populated domains 
form, all of which contribute to the full spectrum of the 
spin excitations, as described below. In the following we 
denote these domains as follows: domain 1 corresponds 
to ki || (1, 1, 1), domain 2 corresponds to k2 || (1, 1, 1), 
domain 3 corresponds to k 3 || (1,1,1), and domain 4 
corresponds to k4 || (1,1,1) (cf. Fig. [I]). We confirmed 
that the domains are equally populated in our measure- 
ments. 

The experimental study of the excitation spectra of the 
helimagnetic state are also extremely challenging because 
of the long wavelength of the helical modulation. This 
makes high-resolution measurements in a relatively small 
portion of reciprocal space around a magnetic satellite re- 
flection defined by q ~ k — 0.035 A^ 1 necessary. In the 



study described here we focus on the cross-over between 
the ferromagnetic limit and the helimagnetic regime de- 
scribed above. In fact, the excitations for even smaller 
momenta deep in the helimagnetic regime and the antifer- 
romagnetic limit described above are technically not ac- 
cessible in present-day inelastic neutron scattering mea- 
surements, even for the most advanced neutron scattering 
techniques. 



II. EXPERIMENTAL METHODS 

A large single crystal of ~ 8 cm 3 grown by the Bridg- 
man method was studied. The paramagnetic spin fluctu- 
ations of the same single crystal were reported previously 
in Ref.^. The single crystal displayed a lattice mosaic of 
the order 0.5°. The specific heat, susceptibility and resis- 
tivity of small pieces taken from this large single crystal 
were in excellent agreement with the literature, where 
the residual resistivity ratio (RRR) was of the order 100. 
The latter indicates good, though not excellent sample 
purity. 

The bulk of our studies was carried out on the cold 
triple-axis spectrometer TASP at the Paul Scherrer Insti- 
tut (PSI) 37 , with a few supplementary measurements on 
PANDA at the Forschungsneutronenquelle Heinz-Maier 
Leibnitz (FRM II). For the experiment the single crystal 
was cooled with an ILL- type orange He-cryostat. The 
sample was oriented with the [1,1,0] and [0,0,1] crystal- 
lographic directions in the scattering plane. In order to 
avoid second order contamination of the neutron beam a 
beryllium filter was inserted between the sample and the 
analyzer. 

The experiments were performed around the nuclear 
Bragg reflection (1,1,0) taking advantage of the large 
magnetic structure factor of the corresponding magnetic 
Bragg reflections (cf. Ref. [33] and Fig. [TJc)). Two dif- 
ferent setups were used. For setup I the triple-axis spec- 
trometer was operated with fixed final wave vector kf 
= 1.2 A -1 . Additionally 20' Soller collimators were in- 
stalled in the incident beam and in front of the analyser 
whereas 40' collimators were used in front of the detec- 
tor. In setup II we used kf = 1.4 A -1 and a collimation 
setting: open-40'-MnSi-40'-80'. The energy resolution in 
setup I and II was approximately 40 /xeV and 90 /ieV, 
respectively. The corresponding Q-rcsolution is depicted 
in Fig. [T] It is important to emphasize that the resolu- 
tion achieved was about the best currently available. 
In order to compare the model with the experimental 
results the theoretical scattering function S(Q,uj) was 
convoluted with the r esolut ion function of the triple-axis 
spectrometers we use d 38 ' 39 -'. For the fits presented in this 
article a program called TASRESFIT^I has been devel- 
oped in PYTHON to perform the convolution of the spec- 
trometer resolution with the theoretical cross-section by 
means of Monte-Carlo integration. 

Further, we have introduced a Gaussian profile centred at 
hoj = 0.01 meV to describe the incoherent scattering to- 
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FIG. 2: Typical constant-Q scans at selected locations of the three trajectories shown in Fig.[T] Note that the energy and 
momentum resolution is at the technical limit currently available. The strong elastic peak at meV is due to incoherent 
scattering. The inset in each panel shows the precise location in reciprocal space as a black spot where data were recorded. The 
curves represent the intensity calculated in the model described in the text, where all data are accounted for by the same values 
of the ferromagnetic spin wave stiffness, the helical wavelength and the intensity, (a), (b): Data for the trajectory with q || ki; 
these scans are dominated by a very broad maximum, (c), (d): Data for q _L ki; these scans show very broad, essentially 
featureless intensity that decreases for increasing energy, (e), (f): Data for an arbitrary trajectory emanating from ki; data for 
this trajectory are characterized by almost featureless intensity over an extremely wide range of energies and a distinct peak in 
a small range. The open data points below -0.4 meV in (f) represent a spurious signal that arises from additional incoherent 
scattering of neutrons from the monochromator crystals of the triple axis spectrometer. 



gether with a constant background for the fit of all data 
presented. The background was carefully measured at 
positions in reciprocal space that were not contaminated 
by inelastic scattering around the (1,1,0) Bragg reflec- 
tion and was estimated as 20 Counts/25 min. 
Finally, some of the energy scans measured during the 
course of our experiments were contaminated by spuri- 
ous scattering that was effectively caused by neutrons 
that were scattered incoherently from the monochroma- 
tor crystals of the used triple-axis spectrometer (see e.g. 
open symbols Fig. |2jf) at E « — 0.4meV). When the 
spectrometer is set to a non-zero energy transfer, i.e. 
|kj| 7^ \kf\, but the angle between k^ and k/ as well 
as the orientation of the sample correspond to a geom- 
etry that allows for Bragg scattering from the sample, 
the incoherently scattered neutrons from the monochro- 
mator that have a wave vector |kji„ c | = |k/| lead to the 
observation of accidental Bragg scattering. This type of 
spurious scattering is well-known as Currat-Axe peaks 
in the literature^} For the analysis we considered only 
scans where the spurious scattering could be separated 
well from the inelastic signal. 



III. EXPERIMENTAL RESULTS 

All scans shown in the following were performed at a 
temperature of 20 K. This way the intensity was strong 
while still being significantly below T c deep inside the 
helical phase. Summarized in Fig.[T|d) are the locations 
in reciprocal space where we have recorded excitation 
spectra. A total of 24 spectra were recorded. To cover 
the regime of interest near the magnetic ordering wave 



vector we have mostly carried out inelastic scans at re- 
ciprocal lattice positions along three traces in reciprocal 
space as shown in Fig. |TJcl) . Two traces (square and 
rhombus symbols) were selected such that they are par- 
allel and perpendicular to fei , while the third trace (circle 
symbols) was recorded along an arbitrary direction. The 
latter may be seen as the most stringent test of the theo- 
retical model used to describe the experimental data. All 
spectra were recorded in energy scans at fixed momen- 
tum Q. The scattering intensity observed experimentally 
is hence the result of contributions from each of the do- 
mains, that is, most of the scattering wave vectors were 
neither perfectly parallel nor perfectly perpendicular to 
the ordering wave vector. 

Prior to our study we expected spectra that are strongly 
reminiscent of the excitations of the ferromagnetic limit. 
In stark contrast, we observed highly anomalous line 
shapes that appear to be inconsistent with any conven- 
tional scenario. For instance, shown in Fig.j2[a) and (b) 
are typical data for q || fc x (squares in Fig.RJ. Here the 
data are characterized by fairly broad dispersive maxima, 
but a naive interpretation of the data suggests an extreme 
form of broadening. This may be compared with data 
recorded for q _L ki shown in Fig.[5|c) and (d) (rhombus 
in Fig.[T]), which is essentially featureless and more like a 
background. Finally, the data for the arbitrary direction 
shown in Fig.|2Je) and (f) display yet another type of 
characteristic, namely distinct maxima suggesting a well 
defined intense mode on a large background of scattering. 
The remarkable variety of seemingly anomalous spectra 
we observe experimentally certainly does not look remi- 
niscent of ferromagnetic spin waves at all. In a first at- 
tempt we tried to fit each scan individually with one or 
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several Gaussians. However, this procedure failed. Rec- 
ognizing that our data was taken for wave vectors much 
smaller than the Fermi momentum, q -C kp, in princi- 
ple, the spectra should be insensitive to most microscopic 
details. We have therefore set up an effective model to 
calculate the excitation spectra based on two physically 
transparent and meaningful parameters, namely the fer- 
romagnetic spin wave stiffness c(20 K ) = 37 meVA 2 de- 
rived from previous experiment J 33 * 35 * and the pitch of the 
helical modulation fch = 0.035 A -1 . A detailed descrip- 
tion of the theoretical model will be given in the follow- 
ing section |IV| To compare the calculated spectra with 
our data we folded them with the energy and momen- 
tum resolution of our measurements (which is about the 
best currently available) and adjusted the absolute value 
of the intensity. We emphasize that the scaling of the 
intensity of all spectra by the same value represents the 
only free parameter of our model. 

Our theoretical model gives complete and precise account 
of all measured features and lineshapes for all momenta 
as illustrated in Fig. [2] for six of the 24 spectra we have 
measured. The agreement between the calculated spec- 
tra and the experimental data can be quantified by the 
remarkably small value of the reduced \ 2 value of ~ 1.3. 
We thereby note that we used a rather conservative es- 
timate of the instrumental resolution, which tentatively 
averages out small details that may possibly exist in the 
data. These underestimates of the experimental resolu- 
tion are quite typical and observed in many inelastic neu- 
tron scattering studies. Further, when artificially mod- 
ifying the values of c and kh by more than ~ 5 % the 
quality of the agreement as measured by \ 2 deteriorates 
rapidly. This suggests strongly that our model sensitively 
captures the entire physics of all excitation spectra mea- 
sured. 



IV. THEORY 

As a starting point of a detailed account of our theoret- 
ical description we emphasize again, that the weakness 
of spin-orbit coupling in MnSi leads to well-separated en- 
ergy and length scales. The underlying physics is based 
on the properties of a fcrromagnet. At energies below 
~ 3 meV, where the spin-flip Stoner continuum sets irP2l, 
and for corresponding small momenta, q < 0.3 A , the 
low-energy excitations of the ferromagnet in the absence 
of spin-orbit interactions are described with high preci- 



(l+Cl+C, (r.l.u.) 

-0.06 -0.03 0.03 0.06 



sion by a simple quadratic dispersion E t 



FM 



cq l . At 



the small wave vectors relevant for our experiment the 
damping of those modes may be neglected. 
Accordingly, the starting point of our theoretical descrip- 
tion is a rotationally invariant non-linear a model (in- 
cluding the appropriate Berry-phase term) to describe 
the low-energy excitations of this ferromagnet. Follow- 
ing Ref. 42 43, we add to this model the leading spin- 
orbit coupling effect, which is given by the rotation- 
ally invariant Dzyaloshinskii-Moriya (DM) interaction 




y 



(H, H, 0) 



(1.042+C,l-042+C, 0.915+0 (r.l.u.) 
-0.12 -0.08 -0.04 0.04 0.08 0.12 

■51 




1x =f 4fc h 




V 



(H, H, 0) 



FIG. 3: Illustration of characteristic features of helimagnon 
bands. In the regime investigated, q ^> fch, a crossover to a 
ferromagnetic dispersion with little weight in side bands was 
generally expected (see text). Instead, for finite q± multi- 
ple bands with approximately equal weight are excited due to 
significant Umklapp scattering, (a) For wave vectors strictly 
parallel to the helix a special symmetry prohibits higher order 
Umklapp processes: a translation of the helix can be com- 
pensated by a rotation of all spins. As a consequence only 
three modes are excited, two modes with minima at ±kh and 
an additional mode centered directly at the position of the 
nuclear Bragg peak with zero momentum but vanishing in- 
tensity for q — > 0. For comparison the dispersion of a fer- 
romagnetic mode to — cq 2 is given (red solid line), (b) The 
strong Umklapp scattering for finite perpendicular momen- 
tum, q± = 4fch, stops the motion of spin excitations with 
small Q|| leading to flat bands (see text and Eq. |6j)). In both 
panels the spectral weight of the corresponding modes is pro- 
portional to the area of the points where we use a maximal 
size for better visibility. For the calculation of the dispersion 
we used c(20 K) = 37 meVA 2 and fc h = 0.035 A" 1 . For clarity 
only a single domain is shown, namely ki . 



Sdm = ff / ^ • (V X $). Being linear in momentum, 
the DM interactions necessarily lead to an instability of 
the ferromagnet and stabilize a chiral helix of the form 



*h(x) = fii cos(k h • x) + n 2 sin(k h • x) 



(1) 



where the unit vectors rii and 112 are perpendicular to 
each other and to the ordering vector kh- In our MnSi 
crystals the sign of g is such that rii, n2 and kh define 
a left-handed coordinate system. Corrections to ([lj are 
small as they arise only from tiny forth order terms in the 
weak spin-orbit coupling neglected in our calculation. In 
a ferromagnetic state these corrections lead to the forma- 
tion of a gap. However, due to the spontaneously broken 
translation invariance and the Goldstone theorem there 
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FIG. 4: Dispersion curves for various directions as calculated 
in the model described in the text. Contributions from differ- 
ent domains are color coded as shown in Fig.[T] The area of 
each point is proportional to the weight of the corresponding 
peak in neutron scattering (only for the most intense peaks 
close to the reciprocal lattice vectors a fixed point size was 
used for better visibility). The complex lineshapes observed 
experimentally (cf. Fig. [2| result from the large number of 
excited bands. The gray shaded bars show the direction and 
parameter range of the constant Q scans recorded in our mea- 
surements, where the width in Q of the gray area represents 
the resolution in Q. Data for the gray bars labeled a through 
d are shown in Fig. [5] (a) Dispersion for the trajectory with 
q || ki. (b) Dispersion for the trajectory with q _L ki. 



cannot be a gap for an incommensurate helimagnet. The 
only effect of the higher-order spin-orbit coupling terms 
is a modification of the spin wave spectrum for momenta 
much smaller than k p <C |kh|, i.e., far b elow the experi- 
mental resolution available to dat d 21 * 22 *. 
To describe the spin-excitations of the helix we expand 
$ around the mean field, <& = $h+5&- As discussed be- 
low the physical situation becomes transparent by using 
a comoving coordinate system, where the order param- 
eter at point x is rotated by the angle kh • x around 
the kh axis using the space-dependent rotation matrix R 
with <&' = In the comoving coordinate system the 
helix is mapped on a ferromagnetic solution <&' h = n\. 
Using length and energy units such that |kh| = 1 and 
-E™ = Ef- h = cfc^ = 1, all free parameters are fixed. 
The long wavelength of the helix with the crystallo- 
graphic lattice (which breaks translation invariance) , im- 



plies a tiny magnetic Brillouin zone and the formation 
of an abundance of bands. For the further discussion it 
is helpful to split the momentum gy parallel to k h into 
q\\ = rtfch + Q\\ introducing a band index n while using 
— < q\\ < fch/2. Denoting by q both qi\ and qj_, 
the Gaussian fluctuations around the mean field are de- 
scribed by the action 

s =2 ^'l q K) M ™'(q,iw Il )(5$Vqk) (2) 

where each entry of M nn i is a two by two matrix describ- 
ing the fluctuations in two spin-directions perpendicular 
to the mean field <&h. For a ferromagnet one needs only 
a single matrix 



M FM (q,c) = 



t 

— iuj q 



2 



(3) 



From the condition detM = one obtains the well-known 
dispersion, oj = q 2 , of a ferromagnet. 
The fluctuations around the helimagnetic state are in 
contrast described by an infinite dimensional matrix, 
where the DM interaction yields two changes with re- 
spect to the ferromagnetic state. First, one obtains an 
additional "1" in the diagonal entries M„ n (q, u>), which 
are given by: 



M nn (q, w) 



-iu) 1 + q 



(4) 



q" = q\ + (<7|| + n) 2 . Second, one obtains non- 
vanishing off-diagonal matrix elements given by: 



M, 



n.nztl 



(q,w) 



-iqf 
iqf 



(5) 



with q^ = q x ±iq y for a helix oriented in z direction. The 



simple matrix defined by Eqns. ( 4|5 | leads to the physics 
explained below and describes the neutron scattering ex- 
periments with high precision and no other parameters 
besides the (measured) pitch of the helix, the measured 
ferromagnetic spin wave stiffness and the overall intensity 
as the only free parameter. 



V. DISCUSSION 

To compare the calculated spectra with the neu- 
tron scattering results, one determines the momentum- 
diagonal 3x3 susceptibility matrix X q (cj) oc (5& (5*) qw 
by inverting first the matrix M and undoing the basis 
change due to the rotation R defined above. The scat- 
tering cross section for neutrons with a given momentum 
transfer Q = G + q close to a reciprocal lattice vec- 
tor G of the MnSi crystalline lattice is proportional to 

(l + n s (w))Im Tr(x q (o; + iO)) - Qxq(^ + *0)Q , where 
ng(w) = l/(e fiu / fcsT ' — 1) is the Bose function and where 
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we have taken into account, that the neutrons couple only 
to the component of \ perpendicular to Q. Within the 
Gaussian approximation one obtains sharp modes. For a 
comparison with experiment these theoretical results are 
convoluted by the experimental resolution. 
We first consider the excitations with vanishing momen- 
tum perpendicular to the propagation vector, = 0, 
predicted in this model. In this limit all off-diagonal 
matrices vanish and one obtains a single mode with the 
dispersion relation ui = cq\\ yjk^ + g 2 , well known from 

the physics of Bose condensed atoms. It descr ibes the 
crossover from a mode with linear dispersio n 1 21 * 32 ! at small 
momentum to a ferromagnetic mode with quadratic dis- 
persion. This crossover is also manifest in the eigen- 
modes as illustrated in Fig.[T] In the ferromagnetic limit 
fc c q» 3> fch the magnetization displays the typical pre- 
cession (Fig. [Tja) , while it oscillates only perpendicular to 
<7H for q {l < fc h (Fig. [TJ;). 

By transforming backto the physical coordinate system, 
one realizes that generically three copies of this mode can 
be observed in neutron scattering with unpolarized neu- 
trons as shown in Fig.|3ja). Besides the two copies with 
minima at ±kjj one also obtains generically a mode di- 
rectly at the Bragg spot with zero momentum. However, 
in contrast to the dispersion of a ferromagnetic mode of 
identical spin stiffness c shown in red, in the helimagnetic 
case the intensity vanishes for q — > 0. The fact that only 
three modes rather than a large number of bands can 
be observed, originates in the special symmetry of the 
helical state: a translation of the helix can be compen- 
sated by an appropriate rotation of the spins. Theoreti- 
cally this property has actually long been knowrJ^l. Only 
higher-order spin-orbit coupling terms which break this 
symmetry may lead to the excitation of further modes, 
however, with a tiny weight. 

The more important second example, which reveals an 
entirely unexpected property, concerns spectra for fi- 
nite perpendicular momentum as shown in Fig. [3j[b). 
Here the off-diagonal terms ^ lead to a mixing of the 
modes. For vanishing qii and q± <C fch one obtainiPIEU 

E{q±) = ^/S/Scq^ which can, however, not be observed 
directly within our experimental resolution. In the limit 
q± fch one finds that more and more modes are ex- 
cited. Mapping the matrix to a harmonic oscillator on a 
lattice (see e.g. Ref. l45|) one finds that typically of the 
order of \J q±/k^ bands are excited with equally spaced 
energies 

E n (q) w c ^ 2 - 2k h q ± + 2yJqJ%(n + ~)\ , 

n = 0,l,2,..V<7±/fch. (6) 

Here two aspects are remarkable. First, for increasing 
perpendicular momentum an increasing number of bands 
are excited. Second, the dispersion is essentially indepen- 
dent of qn for <7|| < \/q±k\ 1 . 

In Fig. plb) we plot the corresponding dispersion as a 
function of q\\ for q± = 4fch with c and fch given above. 



For clarity we consider only the domain associated with 
ki. The transverse momentum strongly couples the fer- 
romagnetic branches (red solid lines), w(q) = c{{q» + 
mfch) 2 + q±) with m = 0, ±1, ±2 . . . , shifted by multiples 
of the reciprocal lattice vector kh of the magnetic Bril- 
louin zone. As the DM interactions grows linearly with q, 
more an more modes are coupled by Umklapp scattering 
when qj_ increases. As described by Eq. Q the repeated 
Umklapp scattering prohibits the motion of the spin ex- 
citations parallel to kh and leads to a very flat dispersion 
for gy < y/q±k h . 

We finally illustrate in Fig. [4] the full complexity of the 
calculated helimagnon bands as probed experimentally 
by the two trajectories in Fig. [I] Here the area of the 
calculated dots is proportional to the intensity and the 
color coding refers to the relevant domain. The gray 
shaded bars show the location and parameter range of 
our measurements, where the width of the gray bars rep- 
resents the resolution of our measurements. Using this 
approach we are able to account for the entire set of very 
broad and extraordinarily rich spectra measured (note 
the very different energy dependences in Fig. [2]) without 
any further adjustments of parameters. 



VI. CONCLUSION 



In many magnetic materials excitation spectra like 
those observed here may be interpreted by invoking a 
host of non-universal properties. For instance, spin- 
orbit coupling and dipolar interactions cause magnetic 
anisotropies and the formation of energy gaps 4 ^. More- 
over, scattering from impurities and crystalline defects, 
the coupling to the particle-hole continuum in metals^, 
magnon-phonorP^ and magnon-magnon interactions^ 
may lead (in combination with spin-orbit interactions) 
to a substantial damping of the magnetic excitations. 
In summary, we report an inelastic neutron scattering 
study of the perhaps simplest manifestation of weak chi- 
ral interactions: the helimagnetic order in the cubic B20 
compound MnSi. The excellent quantitative account we 
achieve establishes the existence of a new universal prop- 
erty of the spin excitations: the formation of a large num- 
ber of dispersive bands of helimagnons. Our study is also 
important for the theoretical description of the skyrmion 
lattice in MnSi which turns out to be stabilized relative 
to a conical phase only by thermal fluctuations. There- 
fore it is important to establish that a quantitative de- 
scription of the fluctuations spectrum can be obtained in 
this material. We expect that similar phenomena will be 
relevant for a broad range of materials where long wave- 
length or incommensurate magnetic order lead to strong 
Umklapp scattering and the formation of multiple bands. 
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